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c§ : ABSTRACT 

H.E.S.S. observed TeV blazar PKS 2155-304 in a strong flare state in 2006 
July. The TeV flux varied on timescale as short as a few minutes, which sets 
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strong constraints on the properties of the emission region. By use of the syn- 
chrotron self-Compton model, we found that models with the bulk Lorentz factor 



Q ■ ~ 100, the size of the emission region ~ 10 15 cm, and magnetic field ~ 0.1 G 



explain the observed spectral energy distribution and the flare timescale ~ a few 
minutes. This model with a large value of T accounts for the emission spectrum 
not only in the TeV band but also in the X-ray band. The major cooling process 
of electrons/positrons in the jet is inverse Compton scattering off synchrotron 
photons. The energy content of the jet is highly dominated by particle kinetic 
energy over magnetic energy. 

Subject headings: BL Lacertae objects: individual (PKS2155-304) - galaxies: 
active - radiation mechanisms: nonthermal 



1. Introduction 

In active galactic nuclei various physical processes such as particle acceleration, emission 
of high energy photons, and so on, take place. Accretion of matter onto the central black 
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holes is thought to be the energy source of these processes. When a fraction of matter is 
ejected from the accretion disks, relativistic jets of plasmas are formed. The relativistic 
jets are thought to explain intense and variable emission from blazars. The beaming effect 
amplifies the radiation from jets depending on the opening angle of the jets, bulk Lorentz 
factor, T, and the angle between the line of sight and the jet axis, 9. Some blazars are 
known to emit very high energy 7-rays i n the TeV ener gy band. They are, for example, PKS 



2155-304, Mrk 421, and Mrk 501 (e.g., IWaener 1120071 ). 



Recent observations of TeV 7-rays have found remarkably short time variability of TeV 
blazars. In 2006 July PKS 2155-304 with redshift z = 0.116 showed an outburst of TeV 
7-rays. The aver age flux during the outb urst was more than 10 times typical values observed 
from the object (jAharonian et al. 1120071 ). During this period X-rays were al so monitored by 



Swift. In the 0.3-10 keV energy band the X-ray flux increased by a factor of 5 (IFoschini et al. 



2007al ). The observations by the High E nergy Stereoscopic Syste m (H.E.S.S.) report that the 
timescale of variation is a few minutes (jAharonian et al. 1 120071 ). The well- resolved burst of 
the TeV 7-ray flux from PKS 2155-304 varied o n timescales only ~ 200 s. Recently MAGIC 
also observed short time variation in Mrk 501 (jAlbert et al. 1120071 ). When the size of the 
emission region in the comoving frame of the jet is denoted by R, the observed timescale of 
the variability sets a limit on R, i.e., R < ct ygx T>/(l + z), where c is the speed of light and 



Aharonian et al. 



(12007( 1 

2GM/c 2 



V= [T(l-p cos 0)]- 1 is the beaming factor with (3 = (1-1/r 2 ) 1 / 2 

argued that V > 60 - 120 R/Rs is required to explain t var ~ 2 min, where Rs 

is the Schwarzschild radius with G and M ~ 1O 9 M being the gravitational constant and 

the central black hole mass, respectively. Such large values of V are also suggested by 

Begelman et al. I (120081 1 recently, based on the requirements of radiative cooling time and 

optical depth for 7-rays due to electron-positron pair production. 

The emission spectra of blazars are characterized by two peaks in the v-vF v representa- 
tion, where F v is the differential energy flux. The lower energy peak is located in the optical - 
X-ray bands and the higher energy peak is located in very high energy 7-ray bands. The lower 
energy peak is most probably by synchrotron radiation of nonthermal electrons/positrons. 
Radiation mechanisms of high energy 7-ray s are thought to be in verse Compton scattering 



off soft photons in the leptonic models (e.g.. lMaraschi et al. 1119921 ) and hadronic interaction 



of re l ativistic particles an d photons in the hadronic models (e.g., iMannheim fc Biermann 
1992c iMiicke et al. 1120031 ). The soft photons of the leptonic models are supplied by syn- 
chrotron radiation by relativistic electrons/positrons in t he jet (synchrotron self- Com pton 



model) or by the external sources such as accretion disks (|Dermer h Schlickeiser Ill9931 ) and 
disk radiation scattered around the jet (ISikora et al. 1119941 ) . These emission models assume 
that the emission region is one zone. 



-3- 



Observations of TeV 7-rays of blazars have revealed that large values of T> > 10 are 
required to fit the emission spectra of very high energy 7-rays by inverse Compton scatter- 

20021 ). On the other hand, the observations of the apparent velocity 



Kino et al. 



ing [e.g., 

of VLBI knots show tha t parsec-scale jets are subrelativistic or at most mildly relativistic 
( jPiner fc Edwards 112004 ). To reconcile the discrepancy between t hese values of T, the decel- 
eration of jets is considered. iGeorganopoulos fc Kazanas I (120031 ) assumed that TeV 7-rays 
are produced by inverse Co mpton scattering off syn chrotron photons that are emitted by 
decelerated jet components. iGhisellini et al. I (120051 ). on the other hand, proposed a spine- 
layer model. In this model they assumed that a fast moving emission region is surrounded 
by a slow moving sheath. Gamma-rays are then produced by inverse Compton scattering off 
the photons emitted in the spine and the sheath. Although various emission models of jets 
with multiple radiation zones have been proposed as mentioned above, a simple synchrotron 
self-Compton (SSC) model is still worth use in obtaining physical parameters of the emission 
regions in jets. 

In this paper we show that a simple SSC model can explain the TeV 7-rays and X-rays 
of PKS 2155-304. PKS 2155-304 is an interesting source because of its strong TeV 7-ray 
emission and the short time variability of TeV emission. Since the redshift of PKS 2155-304 
is 0.1 16, the absorption o f TeV 7-rays by extragalactic background light (EBL) is effective 



Stecker et al. 1120061 ). This is an ideal object to test EBL models as well as emission 



(e.g., 
models. 

In §2] we describe the parameters of our numerical calculations and in £J3] the values of 
the parameters are determined by fitting the observed data. Finally discussion is given in 



Model Parameters 



The emission region is assumed to be a sphere with radius R moving relativistically with 
bulk Lorentz factor T. Below we assume T> = T, i.e., 9 ~ 1/T. We soly e kinetic equations 
of electrons and photons in the emission region (IKusunose et al. Il2000l ). Here we assume 
that the plasmas and radiation in the emission region are in a steady state. Electrons are 
continuously injected in the emission region at rate gWj. The injection spectrum is given by 



9(7) = K-f p exp(— 



7 >7n 



where p, 7 m i n , and 7 max are parameters, and K is the normalization constant determined by 



*nj 



q(j)dj. 



(2) 
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The electrons escape from the emission region by advection on timescale / at j-R/c, where / a d is 
a dimensionless parameter. The cooling processes of the electrons are synchrotron radiation 
and inverse Compton scattering. We assume that the emission region has randomly oriented 
magnetic field B. Photons are assumed to escape from the emission region on timescale 
R/c. We also include the absorption of 7-rays in the emission region by e ± production due 
to photon-photon collisions. This effect is found to be negligible for 7-rays with observed 
energy less than 10 TeV in our numerical results. 

Parameters in our model are T, R, B, g in j, p, 7 m i n , 7max, and / a d- The dependence 
of solutions on / a d is weak and we set / a( j = 4. In addition to the above parameters, the 
cosmological parameters such as Hubble constant H and the density parameters of matter 
Q m and cosmological constant Q\ are needed to calculate the luminosity distance and the 
optical depth for 7-ray absorption by EBL. We assume H = 71 km s^ 1 Mpc -1 , Q m = 0.27, 
and Qa = 0.73. 

It is known that Te V 7-rays emitted by distant sources are absorbed by EBL. This 



was first pointed out by iNikishov I ( 



19621) and d etailed calcula t ions were performed by 



Gould fc Schreder I (119661 ) and I Jelly I (119661 ). Later IStecker et al. I (119921 ) proposed that the 



EBL spectrum is estimated by considering the absorp tion of TeV 7-rays fr om blazars. Since 
then various models of EBL have been proposed (see lStecker et al. 1120061 . for review). The 
redshift of PKS 2155-304 is 0.116 and the optical depth of TeV 7-rays is greater t han unity 



for 7-rays with energy greater than 1 TeV. In our previous work (IKato et al. 1120061 ). we used 



models by lDwek fc Krennrich 



(20 051) to fit the emission spectrum from H1426+428. There 



Dwek fc Krennrich I (120051 ) is applicable in the SSC model. In 



we found that model LLL of 

this paper we use the same EBL model to calculate the deabsorbed TeV spectrum. 



Results 



H.E.S.S. observed a 7-ray outburst from PKS 2155-304 on 2006 July 28 (MJD 53944) 



( lAharonian et al. 1 120071 ) . Almost simultaneous observation by Swift was performed in the 



X-ray band on July 29/30 (IFoschini et al. Il2007at ) . The o bserved emission spectra are shown 
in Figure [TJ The TeV data in 2006 July are from H. E.S.S. (Aharonian et al. 120071 ). The TeV 
spectrum is corrected with an EBL model, LLL, of iDwek fc Krennrich I ( 120051 ) for absorption 
by e production through photon-photon collisions. The X-ray spectrum of 2006 July 29 
an d 30 is shown by a thic k dashed line. This spectrum is a log-parabolic model fit given 
by IFoschini et al. I (J2007bl ). Other X-ray, optical, and radio data are not simultaneous with 
the TeV data. The X-ray data except those of Swift are from BeppoSAX and radio data are 
from NED. 2MASS data are also plotted. The TeV emission spectrum observed in 2003 July 
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( lAharonian et al. Il2005al ) is shown for comparison. 



Our models are shown in Figured] by solid and dashed lines for 2006 July and 2003 July, 
respectively. We did not fit the emission spectra below 10 15 Hz, assuming this emission is 
from different regions, possibly from extended regions far away from the central black hole. 
The values of the parameters for the flare in 2006 July are the following: V = T> = 90, 
B = 0.1 G, R = 9.6 x 10 14 cm, p = 1.9, 7 min = 10, and 7 max = 8 x 10 4 (Model A in Table!]). 
The value of R is ~ 3R$, if M = 10 9 M Q . From these parameters, the timescale of variability 
is t var ~ (1 + z)R/(cF) ~ 400 s. We also obtained the parameter values which give t var ~ 200 
s. Those are shown in Table [T] (Models B and C). When t var ~ 200 s is assumed, the spectral 
energy distribution (SED) of 2006 July is fitted well if F is in the range 100 < T < 150. 
When T < 90, the fluxes in the tails of the lower and higher peaks of the SED are too 
high. On the other hand, when F > 150, the tail of the lower peak (the X-ray band) is too 
steep. SEDs for Models A, B, and C are shown in Fi gure [2] In Table [1 ] the p arameters given 
in other papers such as iFoschini et al. I (2007a ) an d 



Note that t v 



1 h is assumed in 



Foschini et al. I (12007af ) and that iBegelman et al. I (120081 ) 



Begelman et al. I (120081) are also listed . 



did not perform the spectral fitting. 

In Figure [3] various deabsorbed TeV spectra are compared with our models. The deab - 
sorbed spectra are calculated with different EBL models given by lDwek fc Krennrich I (120051 ). 
When EBL models other than LLL, LHH, LLH, and LHL are applied (LLH and LHL are 
not shown in the figure), the TeV spectrum has a peak at ~ 3 x 10 26 Hz. Such a spectrum 
is difficult to produce by the one-zone SSC model, if the X-rays are emitted in the same 
region as TeV 7-rays. In particular the Klein-Nishina effect suppresses the emission in the 
TeV band. 

In our simulations, small numbers of nonthermal electrons and synchrotron photons are 
initially injected and time evolution is followed until a steady state is attained. Because of 
this initial condition, the SED takes longer time than the observed timescale to attain a 
steady state. The time evolution of the SED of Model A is shown in Figure H] In actual 
flares, the initial condition may different from that used here and flares are not in a steady 
state. We used a steady state to obtain typical values of the source parameters. 

In Figure [5] the electron spectrum for the flare is shown. The electron spectrum has a 
cutoff at 7 ~ 10 4 . Because the injection spectrum has an exponential cutoff at 7 max > 10 4 , 
efficient Compton cooling made the cutoff energy smaller. The synchrotron cooling time of 
electrons with 7 = 10 4 is about 7.7 x 10 6 s for B = 0.1 G. The energy density of synchrotron 
photons is w syn ~ 5.3 x 10~ 3 ergs cm -3 in Model A. This results in the Compton cooling 
time ~ 5.8 x 10 5 s for electrons with 7 = 10 4 . Electrons with 7 ~ 7max are cooled in the 
flare timescale in the comoving frame of the jet. The Compton cooling time is short enough 
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to be a major cooling process during the TeV flare. 

The numerical results show that the energy densities of nonthermal electrons and mag- 
netic fields in Model A are, respectively, u e = 3.8 ergs cm -3 and ub = 4.0 x 10~ 4 ergs 
cm -3 . This yields a large value of u e /uB ~ 9.5 x 10 3 . The values of energy contents of 
different models are shown in Table [2j The dominance of particle kinetic energy over mag- 
netic energy should be taken into account in consideration of acceleration mechanisms. The 
powers contained in the jet are, in the form of radiation and electron kinetic energy, given 
by L rad = iicR 2 T 2 u Ta d and L kin = iicR 2 T 2 u e , respectively, where w ra d is the radiation energy 
density. The numerical results are given in Table [2j 

Our model include internal absorption of 7-rays by e^ production. However, our nu- 
merical results show that the 7-ray absorption is negligible. Because of the large values of 
T, the soft photon density and the maximum 7-ray energy in the comoving frame of the jet 
are not large enough for e ± pair production. 

In FigureHJwe show a TeV spectrum observed in 2003 July as well. In 2003, TeV emission 



was observed by H.E.S.S. several times (jAharonian et al. Il2005al ). In those observations the 



TeV flux was much lower than that in the 2006 July flare. The TeV spectrum of 2003 July 
are fitted with parameters such as T = T> = 40, B = 0.08 G, R = 2.5 x 10 16 cm, p = 1.95, 
7min = 10, and 7 max = 2 x 10 5 . With these parameters, numerical results give u e /uB = 26.7. 
Compared with the flare of 2006 July, a blob with a larger size but a smaller value of V was 
involved. The powers of radiation and electrons are, respectively, L ra d = 1.7 x 10 43 ergs s" 1 
and Lkin = 6.4 x 10 44 ergs s~ 4 . Compared with the flare state in 2006, L^ n is smaller by a 
factor 0.24 and L ra d is smaller by a factor 0.52, if Model A is adopted. 



In lAharonian et al. I (j2005bl ). they derived source parameters for the H.E.S.S. obser- 
vation in 2003, October and November, with different EBL models from that used here. 
According to their results, the source was in a low or quiet state. With t heir leptonic mode l, 



they obtained R = 1.5 x 10 15 cm, V = 25, and B = 0.25 G (Model 2). iKino et al. I (120021 ). 
on the other hand, used the SSC model to fit another quiescent state, but without the cor- 
rection for 7-ray absorption by EBL. They obtained R = 9 x 10 15 cm, T> = 33, and B = 0.3 
G. The ratio u s /ub is 3, which is smaller by a factor 9, compared with our results for 2003 
July. The source size obtained by us is the largest among these models for quiescent states. 



4. Discussion 

In this work we have shown that a large value of the bulk Lorentz factor, T ~ 100, is 
required to explain the emission from PKS 2155-304 in the TeV band as well as the X-ray 



-7- 



band. With R = 9.6 x 10 14 cm and T = 90, the timescale t var ~ 400 s is obtained, while 
t var ~ 200 s is obtained for R ~ 5 x 10 14 cm and r ~ 100 - 150. The value of V = T 
as large as 100 requires the viewing angle as small as 0.01 rad, which makes the case an 
unlikely coincidence if the opening angle of the jet is the same order as 1/T. However, if the 
opening angle is larger than 1/T, we can observe a portion of the jet that aligns with the 
line of sight, thus avoiding an excessively small observability. Considering that Mrk 501 also 
exhibited a similar short timescale TeV flare, we think that the opening angle is larger than 
1/T (by a factor a). If this is the case, R represents the lateral size along the jet motion, 
while the transverse size is larger than R by the same factor a. Thus the real size of the 
emission region is larger than R and only a part of the emission regions is observed owing to 
a strong relativistic beaming effect. Then the kinetic power should be by a factor a 2 larger 
than a conventional estimate given by ircR 2 T 2 u e . A further speculative possibility is that 
the opening angle is as small as 1/T but the jet axis wanders within a large angle and the 
flare is observed when its direction happens to fall in the line of sight. 

Our results show that the jet plasma is highly particle energy dominated over magnetic 
energy, i.e., u p /u r ~ 10 4 in the flare state in 2006. The ratio u e /uB has been discussed in 



many papers. In iKino et al. I (120021 ) we discussed that TeV blazars are generally particle 
dominated typically by a factor of 10. For a quiescent state of PKS 2155-304 we obtained 
the ratio of 3. Thus the ratio of 10 4 seems unusually high, although there is no serious 
theoretical reason to expect a ratio of unity except that it leads to a minimum kinetic 
power for a given observed spectrum. The obtained high value of the ratio basically results 
from the fact that the SSC luminosity is much larger than the synchrotron luminosity. In 
this respect, we note that the X-ray and TeV observations are not strictly simultaneous; 
if the X-ray flux during the TeV flare is higher than the adopted value in this paper, the 
ratio decreases correspondingly. Also, some soft photon sources other than synchrotron 
radiation will help to decrease the ratio. However, considering a good fit to the observed 
spectrum obtaine d here, such poss ibilit ies are not easy t o rea lize and the reduction will 



not be significant. iLevinson I (120071 ) and lBegelman et al. I (120081 ) adopted the view that the 
TeV flare is caused by external Compton mechanism. In this case, the model may avoid 
an extremely large value of u e /uB- However, considering the short timescale of radiative 
deceleration, the acceleration of the emitting plasma to the large bulk Lorentz factor itself 
becomes very difficult because of the large radiation drag effect. Also note that both papers 
did not make detailed multi-wavelength fitting. 



As mentioned above, iBegelman et al. I (120081 ) argued that inverse Compton scattering 
off external soft photons is favored as the dominant radiation process to produce TeV 7-rays. 
The major differences between their model and ours are that they assumed the dominance 
of the Poynting flux over the kinetic energy flux as well as the synchrotron radiation as a 



major cooling process to estimate particle's cooling timescale, while in our model the energy 
flux is dominated by nonthermal particle's kinetic energy and the cooling is dominated by 

SSC. 

In comparison between the flare state in 2006 July and the quiescent state in 2003 July, 
the powers in radiation and electrons increased by factors 2 and 4.2, respectively, during the 
flare in 2006, if Model A is adopted. The increase in the jet power is significant and the 
efficient acceleration of electrons is implied. The size of the blob for the 2006 July flare is 
~ 10 15 cm. This is about 3 R$, if the mass of the central black hole is 1O 9 M . A large 
energy injection into a very compact region occurred in the flare. The obtained parameter 
values for the flare event and various phases are not well ordered and we found no simple 
tendencies from quiescent states to flares. The TeV flare may be caused by a concentration 
of energy into a sub-horizon scale region, while during quiescent states we may be seeing 
larger scale energy dissipation on longer timescales. 



In applying our model, we assumed a low flux model of EBL, i.e., LLL of lDwek fc Krennrich 



(120051 ) . With this EBL model the one- zone SS C model is found to be via b le for this 



TeV blazar, although other models proposed by iGeorganopoulos fc Kazanas I (120031 ) and 



Ghisellini et al. I (120051 ) might explain the TeV spectrum with higher fluxes of EBL. This 



should be studied in future work. 

Radiation below 10 15 Hz is not fitted with our model. We assumed that the emis- 
sion comes from more e xtended regions far away from the central black hole (see also 



Aharonian et al. I l2005bl ). If the radiation below 10 15 Hz is emitted by the jet plasma, 
inverse Compton scattering produces radiation in the MeV - GeV band. Observation in 
this energy band is particularly important to set constraints on the emission mechanisms of 
blazars. 
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from the Ministry of Education, Culture, Sports, Science and Technology of Japan. 
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Fig. 1. — The emission spectra of PKS 2155-304. The TeV data are observed in 2006 
July (MJD53944) and 2003 July. The data of MJD53944 (open squares) and those of 2003 
July are corrected for absorption by EBL. The thick dashed line is X-ray data for 2006 July 
observed by Swift. Other X-ray data are obtained by BeppoSAX, which are not simultaneous 
with the TeV data. Model A is shown by a solid line and a model for 2003 is shown by a 
dashed line. 
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Fig. 2. — SEDs of Models A, B, and C are shown by solid, dashed, and dot-dashed lines, 
respectively. 
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Fig. 3. — The TeV spectrum of 2 006 July (filled circles) is deabsorbed with different EBL 
models (IDwek fc Krennrich 1 120051 ) . Models A, B, and C are shown by solid, dashed, and 
dot-dashed lines, respectively. 
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Fig. 4. — Time evolution of SED of Model A is shown. The SEDs are plotted at every R/c 
and evolve from lower to upper lines. Labels are t/(R/c) in the comoving frame. 
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Fig. 5. — The electron spectrum for the flare in 2006 July. Models A, B, and C are shown 
by solid, dashed, and dot-dashed lines, respectively. 
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Table 1. Model Parameters 



Parameter Model A Model B Model C Foschini et al. Model 8, Begelman et al. Model b 

>50 

> 1.6 

4.5 



r 


90 


100 


120 


30 


V 


90 


100 


120 


33.5 


B(G) 


0.1 


0.14 


0.08 


0.27 


R (10 14 cm) 


9.6 


5.4 


6.5 


50 


V 


1.9 


1.9 


1.9 


2.5 C 


7max (10 4 ) 


8.0 


6.7 


5.5 


17.5 d 


7min 


10 


10 


10 




/break 








1.5 x 1C 


/peak 











10 4 



a The model for July 29 given in iFoschini et al. I (J2007af ). The variability timescale ~1 h was 
assumed. 



b The model given in iBegelman et al. I (120081 ) . T = 50, t var = 300 s, and isotropic luminosity 
L iso = 10 46 ergs s _1 are assumed to calculate B, R, and 7 pea k, where 7 pea k is the electron Lorentz 
factor that emits synchrotron radiation peaking at ~ 10 16 /r Hz. 

c The power law index for 7 min < 7 < 7 br eak- The index for 7 bre ak < 7 < 7max is p + 1 

d The maximum energy of electrons. 

e The break energy of electrons obeying a broken power law. 
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Table 2. Energy Contents 



Parameter 


Model A 


Model B 


Model C 


r 


90 


100 


120 


u e (ergs cm -3 ) 


3.8 


8.7 


5.7 


Ub (1CT 4 ergs cm -3 ) 


4.0 


7.8 


2.6 


u e /u B (10 4 ) 


0.95 


1.1 


2.2 


L rad (10 43 ergs s" 1 ) 


3.3 


2.9 


1.7 


L kin (10 45 ergs s" 1 ) 


2.7 


2.4 


3.3 



